Scaling of cosmic string loops 
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We study the spectrum of loops as a part of a complete network of cosmic strings in flat spacetime. 
After a long transient regime, characterized by production of small loops at the scale of the initial 
conditions, it appears that a true scaling regime takes over. In this final regime the characteristic 
length of loops scales as O.lt, in contrast to earlier simulations which found tiny loops. We expect 
the expanding-universe behavior to be qualitatively similar. The large loop sizes have important 
cosmological implications. In particular, the nucleosynthesis bound becomes G/i < lO"'^, much 
tighter than before. 
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Cosmic strings could be formed as linear defects at 
symmetry breaking phase transitions in the early uni- 
verse 0. Alternatively, they could arise as fundamen- 
tal or j6-strings at the end of brane inflation Qi d, j^. 
Strings can produce a variety of observational effects, and 
searches are now underway for their signatures in gravita- 
tional lensing, CMB anisotropics, and gravitational wave 
background. A good theoretical understanding of string 
networks is crucial for interpreting the results of these 
searches. However, despite much effort, the evolution of 
cosmic strings is not yet fully understood. 

An evolving string network consists of two components: 
long strings and sub-horizon closed loops. It is fairly 
well established that long strings exhibit scaling behav- 
ior: both the average distance between the strings d{t) 
and the coherence (or persistence) length ^(t) scale with 
the cosmic horizon. 



d{t) ^ m ~ t. 



(1) 



In the early work on cosmic strings, it was expected 
0] that the typical length of closed loops scales in a 
similar manner, l{t) ~ t. The first numerical simulations 
of string evolution seemed to support this scenario 
However, later, high-resolution simulations 0,0 showed 
that the loop sizes were actually much smaller than the 
horizon and gave no evidence for scaling. On the con- 
trary, the typical loop size remained nearly constant and 
close to the resolution limit of the simulations. The sim- 
ulations also revealed that long strings had a significant 
substructure, with short-wavelength wiggles all the way 
down to the resolution limit. 

The standard scenario of string evolution that has 
emerged from these findings assumes that the typical size 
of loops l{t) is set by the scale of the smallest wiggles, 
which is in turn determined by damping due to gravita- 
tional radiation The typical loop size is then given 
by 



l{t) - at 



(2) 



with a ~ {G^jl)^ , where G is Newton's constant and 
/i is the mass per unit length of string. With plausible 



assumptions about the spectrum of wiggles, the power in- 
dex 7 is in the range 2 < 7 < 3. The observational bound 
on the string mass parameter Gji is G/x < 10~^, and the 
corresponding bound on a is a < 10~^^, indicating that 
the loops are extremely small. 

A more radical string evolution scenario, proposed in 
[l^ . suggests that the loops are even smaller. It claims 
that strings lose most of their energy by direct emission of 
microscopic loops of size not much greater than the string 
thickness. This idea was not confirmed in other simula- 
tions 0, It is hard to tell which, if any, of these 
scenarios is correct, since the loop sizes they suggest are 
well beyond the resolution limits of current simulations. 

To address this issue, we developed in a flat-space 
string simulation which does not suffer from the problem 
of smallest resolution scale. This simulation uses func- 
tional forms for the string positions and is exact to the 
limits of computer arithmetics. In Ref. |lf)l | we used our 
simulation to show that the spectrum of wiggles on long 
strings scales with time, even in the absence of gravita- 
tional damping. The spectrum has a universal form; it 
is peaked at / ^ 0.3i and declines slowly towards smaller 
scales. 

In the present paper we report on the study of the 
closed loop component of the simulation. We find that 
loops are produced in a wide range of sizes, with most of 
the string length going into relatively large loops of size 
I ^ O.lt, comparable to the inter-string distance. This 
scaling behavior is established only after a long transient 
regime, characterized by copious production of tiny loops 
whose size is set by the scale of the initial string network. 
We believe that it was this transient regime that was 
observed in earlier simulations. 

The numerical results presented here are based on the 
computer simulation described in . In order to be able 
to study the evolution at late times we periodically in- 
crease the simulation volume as described in |0 . Specif- 
ically, we create 8 identical replicas of the simulation box 
and glue them together to form a box twice the size of 
the original one. The resulting unphysical correlations on 
super-horizon scales are removed by allowing the string 
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reconnection probability to be p < 1. We used the value 
of p = 0.5 which maximizes the rate of decay of the corre- 
lations.^ We have verified that simulations with p = 0.2 
and 0.8 give similar results. 

Until recently, there was a widespread opinion that all 
values of p other than p = \ are unphysical and there- 
fore of little interest. However, it was pointed out in 
[^L Ufillr^ that for fundamental and ZJ-stringsp generally 
differs from 1, and may even be <^ 1. The dependence 
of network evolution on the value of p is an important 
problem, but we will not try to address it here. It will 
require more extensive simulations, since the evolution is 
very slow for small p. 

To start our simulation we generate four different 
Vachaspati-Vilenkin initial conditions ^3 and overlay 
them in the same simulation box. We displace the four 
realizations relative to each other so that equivalent lat- 
tice points lie on the corners of a tetrahedron of height 
0.5. In previous simulation work we found that the corre- 
lation length ^{t) is larger than the inter-string distance 
d{t) by a factor of about 2.6 in the scaling regime, while 
in the initial conditions this factor is only about 1.3. By 
overlaying four initial conditions, we decrease the initial 
inter-string distance by a factor 2, while not changing the 
correlation length, so the initial conditions have proper- 
ties more similar to a scaling regime. 

The maximum initial box size with four interlaced ini- 
tial conditions that we can simulate on a single processor 
with 3GB of memory is 50, which becomes 800 by time 
1000, after 4 doublings of the box size. The results pre- 
sented here are the averages of 25 simulations. The ex- 
pansion of the simulation volumes takes place at typical 
times 33, 93, 259, and 924. 

A snapshot of a cubic section of the network at time 
800 can be seen in Fig. ^ There are some long strings 
crossing the cube, many small loops, and a few loops 
of sizes up to 60, comparable to the average inter-string 
separation at that time. 

The spectrum of wiggles, as defined in 0, is shown 
in Fig.|21as a function of kt. There are two peaks present 
on the plot: one scaling and one non-scaling. The non- 
scaling peak slowly decays and moves to large values of 
fct, but stays roughly at the same value of k that cor- 
responds to the initial correlation length. In contrast, 
the scaling peak remains at almost the same value of kt. 
Note that in the non-scaling peak was completely 
eliminated by smoothing. We have not used smoothing 
here, because it potentially distorts the spectrum of small 



To minimize disturbance to the string network we do not use 
the smoothing procedure of Il5l . and we keep the reconnection 
probability always 0.5, rather than switching between 0.5 and 1. 
Apart from making the time dependence smoother, this has no 
significant effect on the results. 




FIG. 1: A 150^ section of the network at time 800. Loops of 
sizes less than 10 are not shown. 
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FIG. 2: Evolution of the power spectrum P{kt) vs. kt at 
times: 60, 100, 160, 250, 400, 600, 900 (bold). Error bars are 
the one sigma run-to-run variation. 



loops and because the jumps it introduces make it diffi- 
cult to see trends in the loop production spectrum. 

We characterize the rate of loop production by the 
function n(Z, t) — the number of loops produced per unit 
loop length per unit volume of the network per unit time. 
In a scaling network, the number of loops with sizes be- 
tween I and I -\- dl produced in a volume evolving from 
time t to t -\- dt is the same as the number with sizes be- 
tween 21 and 21 + 2dl produced in a volume 8L'^ evolving 
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FIG. 3: The primary loops production function fp{x) in 
equal logarithmic bins of time centered at times 64, 89, 125, 
175, 245, 342, 479, 671, and 939 (bold). 

from time 2t to 2t + 2dt. Thus, for scaling, 

n{l,t)=t-''f{x), (3) 

where / can be any function of a; = l/t. 

In a cosmological string network, infinite strings self- 
intersect and produce loops. These loops can ether 
reconnect with other strings, fragment further by self- 
intersections, or oscillate without self-intersections. But 
in a simulation, all strings are closed, so we need some 
definition of the point at which a loop is considered to 
have been produced. We proceed as follows. 

First, we say that a loop is a survivor if neither it, 
nor any fragment produced from it, rejoins any other 
string. To conserve computer memory, we do not allow 
strings shorter than a minimum length nt to rejoin the 
network, so any string smaller than nt is automatically 
a survivor. We have verified in [isj that the network 
evolution is rather insensitive to the choice of k, and is 
not significantly modified even if one sets k = 0. In the 
present simulation we used k, = 0.25. 

A loop is primary if it is a survivor but none of its an- 
cestors are survivors. In Fig. |21we plot the primary loop 
production function x^fp{x), where fp{x) = t^np{l,t) 
and Up is the production function of primary loops with 
the same conventions as for n. The graph thus shows 
the fraction of total length produced in primary loops 
for each logarithmic interval in x. 

For a scaling spectrum, we expect fp{x) to be indepen- 
dent of time. Instead, we see two peaks closely related to 
the two peaks in the power spectrum. The scaling peak 
(on the right) does not move in x but remains at 

Ip - OM (4) 

and increases in amplitude. On the other hand, the non- 
scaling peak (on the left) always remains at the scale 
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FIG. 4: The final production function of loops x^f{x) for the 
same time intervals as Fig. |21 

of the initial correlation length, so it moves to smaller 
values of x, and decreases in amplitude. At early times 
the production of loops is dominated by short lengths, 
but as long strings become smoother, the loop production 
at small scales decreases. In contrast, the scaling peak is 
sub-dominant in the beginning, but as the time advances 
it steadily grows. At even later times, we expect the peak 
related to initial conditions to vanish, leaving a scaling 
spectrum. 

A primary loop of size I is likely to fragment in time 
t ~ 1/4. The fragmentation process continues until all 
loops find themselves in non-self-intersecting trajectories. 
In Fig. 0] we plot the final production function of loops 
x'^ f{x). As for primary loops, there are two peaks in the 
final loop production function: one scaling and one non- 
scaling. The non-scaling peak is going down and we ex- 
pect it eventually to vanish, while the scaling peak slowly 
becomes the dominant one and remains at constant x, 

If^G.lt. (5) 

The non-scaling peak in Fig. ^ is present only because 
of the production of small primary loops shown in Fig. |31 
If one considers only final loops produced from primary 
loops whose length is greater than 8.4, the small-scale 
peak is absent, as shown in Fig. El = 8.4 is about 
twice the size of the smallest loops in the initial string 
network. This particular choice is due to the logarithmic 
binning we used for the loop statistics.) 

The scaling part of the final loop production function 
can be fit by a power law 

/(x)«Ax-^ (6) 

with 

/? = 1.63 ±0.03. (7) 
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FIG. 5: The final production function of loops x'^f{x) from 
the last time interval in Fig.0] considering only primary loops 
longer than 8.4. 



We found A = 82 ±2 in our simulations, but since a large 
part of the loop production is still at small sizes even 
at the last times simulated, we expect the final scaling 
value of A to be somewhat larger. (The error bars above 
indicate only statistical variation between different runs 
of the simulation.) 

The corresponding loop production rate is: 



n{l,t) ^ t-'^ fix) ^ At/^^^l-f^ , 



(8) 



with a sharp cut-off for large loops at I ^ O.lt. From 
Eq. the total number of loops produced is divergent, 
while their total length is finite. 

Our simulations are done only in flat space, but one 
can make a reasonable conjecture about the expanding 
universe. The expansion of the universe stretches the 
excitations on the strings, so small-scale structure tends 
to be reduced relative to the flat-space case. Thus we 
expect the non-scaling peak at the initial condition size 
to be eliminated more quickly, and the loop production 
spectrum to fall more rapidly toward small sizes in the 
expanding universe than in flat space. If as a result the 
power index (3 is decreased by at least 0.63, Eq. (jSJ will 
give a convergent total number of small loops. 

In any case, we expect, on the basis of the simulations 
described here, that most of the energy produced in loops 
from a cosmic string network appears at scales compa- 
rable to the inter-string distance. This agrees with early 
expectations [^|^, but not with later simulations 



If there is a scaling process of loop production in an 
expanding universe, then loops produced over time at 
some fixed fraction a of the horizon size give rise to a 
spectrum of presently existing loops 



N{l,t) oc l-/^" 



(9) 



where f3c = 5/2 for a radiation-dominated universe and 
Pc — "2 for matter. This spectrum rises more steeply to- 
ward smaller I than the loop production function, Eq. 
©, indicating that late loop production does not signif- 
icantly affect the form of the spectrum, Eq. Thus as 
long as we know that /? < /?c in the expanding universe, 
we do not need to know the precise value of /3. The form 
of the loop distribution is given by Eq. 0. This dis- 
tribution is cut off at small scales by gravitational back 
reaction. 

The large value of a suggested by our simulations im- 
plies a higher energy density of the string-generated grav- 
itational wave background and tighter observational con- 
straints on the string parameter Gfi. The requirement 
that gravitational waves from strings do not affect the 
predictions of big bang nucleosynthesis can be expressed 
as liej (for a > lO^Gn) 



(10) 



It follows from the analysis in [23 that for a ^ lO^G/i 
the millisecond pulsar observations yield the bound Gfi ^ 
10^'^, which is similar to ifTIHl for a 0.1. (Both of 
these bounds assume that the reconnection probability 
is p ~ 1. We expect p to produce a denser string net- 
work and so make this bound more stringent. Unfortu- 
nately simulation of small p is difficult because the net- 
work evolves much more slowly.) 

Recent work by Ringeval, Sakellariadou and 

Bouchet (RSB) discussed loop production in expanding 
universe simulations. They found that the distribution of 
loops grows steeply toward small scales, with the power 
index /? « 3.0 in the radiation era and /? ~ 2.5 in the 
matter era. Since these indices are larger than (3c, the 
index /3 in their loop production function must also have 
these values, in sharp contrast with our result, Eq. 0. 
The most likely explanation of this discrepancy is that 
the RSB simulation is still in the transient regime dom- 
inated by small loop production. Indeed, the shape of 
the loop production function obtained by RSB is similar 
to that in our Fig. 0| at times t < 100, when the scaling 
peak is not yet pronounced. 

This work was supported in part by the National Sci- 
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